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Abstract. Solutions to the wave equation on de Sitter-Schwarzschild space 
with smooth initial data on a Cauchy surface are shown to decay exponentially 
to a constant at temporal infinity, with corresponding uniform decay on the 
appropriately compactified space. 



1. Introduction 

In this paper we describe the asymptotics of solutions to the wave equation on 
dc Sitter-Schwarzschild space. The static model for the latter is M = Mf x X, 
X = {rhhiTds)r X with the Lorentzian metric 

(1.1) g = fidt'^ - fi'^ dr^ -r^ du^, 
where 

(1.2) — 

with A and m suitable positive constants, < 9m^A < 1, rbh, ''ds the two positive 
roots of /i and duj"^ the standard metric on . We also consider the compactification 

of a: to 

X = [rbh,rds]r X S^. 
Then /Lt is a defining function for dX since it vanishes simply at rbh,7'dS: i-e. 2/3 = 
^ ^ at r = rbh,''dS- Moreover, in what follows we will sometimes consider 

(1.3) a = ^5 

as a boundary defining function for a different compactification of X. This amounts 
to changing the C°° structure of X by adjoining a as a smooth function. Wc denote 
the new manifold by Xi . 

The d'Alembertian with respect to (|l.ip is 

(1.4) □ = a~^{D'^ - a'^r-^Drir^a'^Dr) - a^r^^A^,), 

where A^^ is the Laplacian on S^. We shall consider solutions to Du = on M. 

Regarding space-time as a product, up to the conformal factor a^, is in fact 
misleading in several ways - in particular, solutions to the wave equation do not have 
simple asymptotic behavior on this space. Starting from the stationary description 
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of the metric, it is natural to first compactify the time hne exponentially to an 
interval [0, 1]t- This can be done using a diffeomorphism T : M — > (0, 1) with 
derivative T' < 0. Set 

(1.5) T+ = Ta,+ = mt>C, 
with A to be determined and let 

(1.6) T = T+mt>C. 
Similarly set 

(1.7) T_ = Ta,^ = e^^*, T = 1 - r_ in t < -C. 

Near infinity T depends on the free parameter A. The boundary hypersurface T+ = 
(i.e. T = 0) in 

(1.8) [0,1]txX 

is called here the future temporal face, T_ = the past temporal face, while r = rbh 
and r = r^s arc the black hole, rcsp. de Sitter, infinity, or together spatial infinity. 

In fact, it turns out that we need to use different values of A at the two ends, Abh 
and Ads- This is discussed in more detail in the next section. There are product 
decompositions near these boundaries 

[0, 1]t X [rbh,rbh + ^) x S^, [0, 1]t x (ras - S,rds] x 

If S is so large that those overlap, the transition function is not smooth but rather 
is given by taking positive powers of the defining function of the future temporal 
face, so the resulting space should really be thought of having a polyhomogeneous 
conormal (but not smooth) structure in the sense of differentiability up to the 
temporal faces. In particular, there is no globally preferred boundary defining 
function for the temporal face, rather such a function is only determined up to 
positive powers and multiplication by positive factors. Thus, there is no fully 
natural 'unit' of decay but we consider powers of e~*, resp. e*, in a neighborhood 
of the future and past temporal faces, respectively. 

It turns out that there are two resolutions of this compactificd space which play 
a useful role in describing asymptotics. The first arises by blowing up the corners 

{0} X {rbh} X {0} X {rds} x §^ {1} x {rth} x {1} x {rds} x S^^ 

where the blow-up is understood to be the standard spherical blow-up when locally 
the future temporal face is defined by Tajs,+ , resp. TA^h,-!- at the de Sitter and 
black hole ends. The resulting space is denoted M . The lift of the temporal and 
spatial faces retain their names, while the new front faces are called the scattering 
faces. This is closely related to to the Penrose compactification, where however the 
temporal faces arc compressed. 

Thus, a neighborhood of the lift of {0} x {rbh} x §^ is diffeomorphic to 

(1-9) [0,e)p X [rbhjJ-bh + (5) X p = pbh.+ = T^Ath.+Z^- 

Similarly, a neighborhood of {0} x {rds} x §^ is diffeomorphic to 

(1.10) [0, e)p X (rds - <5, rds] x S^, p = pds,+ = ^a.^.+Zm- 

If (5 > is large enough, these cover a neighborhood of the future temporal face 
tf+, given by the lift of T = 0. Thus a neighborhood of the interior of tf+, is 
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polyhomogeneous-difFcomorphic to an open subset of 
[0,e)^ X (rbh,r-ds) x §1, 

^^■^^^ l/(2Abh. + ) r l/(2AdS +) f 

^ = Phh,+ foi' near r^i, x = p^g^^ for r near r^s, 

where we let the preferred defining function (up to taking positive multiples) of tf + 
be a; = e^* in the interior of tf+, hence x = pj,^*'^^'''''^'' at the black hole boundary 
of tf+. This means, in particular, that a neighborhood of tf+ is polyhomogeneous 
diffeomorphic to 

(1.12) [0,e), X [rbh,rds] X 

If fj, is replaced by a as the defining function of the boundary of X, i.e. X and 
r+ are replaced by X1/2 and T_^' (and analogously in the past) the resulting space 
is denoted Mi/2. Thus M1/2 is the square-root blow up of M, where the square 
root of the defining function of every boundary hypersurface has been appended to 
the smooth structure. Here tf+ is naturally diffeomorphic to X in M, and to X1/2 
in Afi/2- Both M and M1/2 have polyhomogeneous conormal structures at tf+ and 
tf_; we let the preferred defining function (up to taking positive multiples) of tf+ 
be a; = e~* in the interior of tf+, hence x = p^l^"^^^ at 9tf+. 




Figure 1. On the left, the space-time product compactification of 
de Sitter-Schwarzschild space is shown (ignoring the product with 
S^), with the time and space coordinate lines indicated by thin 
lines. On the right, M is shown, with the time and space coordi- 
nates indicated by thin lines. These are no longer valid coordinates 
on M . Valid coordinates near the top left corner are p and /i. 

Solutions to the wave equation, when lifted to this space have simpler asymp- 
totics than on the product compactification, (jl.Sp . The first indication of this is 
that g extends to be C°° and non-degenerate, up to the scattering faces, p = 0, 
away from spatial infinity and uniformly up to the temporal face; the scattering 
faces are characteristic with respect to the metric. We can thus extend M across 
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/i = to a manifold M, by allowing jj, to take negative values; then the scattering 
face becomes an interior characteristic hypersurface. 

A further indication of the utility of this space can be seen from our main result 
which is stated in terms of 

(1.13) ^;?+(M). 

This consists of those functions which arc on M away from tf_(-, and are conormal 
at tf+, including smoothness up to the boundary of tf_(-. Such spaces are well- 
defined, even though the smooth structure on M is not; the conormal structure 
suffices. Thus, the elements of (11.131) . are fixed by the condition that for any k and 
smooth vector fields Vi, . . . ,Vk on M which are tangent to tf+, 

Vi.-.Vkvex'^Ll.t^iM), 

where (M) is the L^-space with respect a density z^b such that xi^h is smooth 

and strictly positive on M. Such a density is well-defined up to a strictly positive 
polyhomogeneous multiple even under the operation of replacing a; by a positive 
power, although the weight is not. Thus, for all to g R, 

x"+'C°"(M) c A'S^iM) c x"'L°°{M), e > 0. 

The main result on wave propagation is: 

Theorem 1.1. Suppose u e C°°{M) satisfies Du = for x £ (0,1), then there 
exists a constant c and e > such that 

u~cGAt,^{M)^x'A%^{M). 

Thus, u has an asymptotic limit, which happens to be a constant, at tf_|_, uni- 
formly on X. 

While we have concerned ourselves with the behavior of the metric at the corner, 
in regions where p < C (i.e. near temporal infinity), it is worthwhile considering 
what happens where p > C , i.e. at spatial infinity. As we shall see, spatial infinity 
can be blown down, i.e. there is a manifold M and a map /3, /3 : M ^ M such 
that /3 is a diffeomorphism away from spatial infinity, and such that g lifts to a C°° 
Lorentz b-metric on M, with tangent (i.e. b-) behavior at the temporal face, smooth 
at the other faces, with respect to which the non-temporal faces are characteristic. 
One valid coordinate system in a neighborhood of the image of a neighborhood of 
the black hole end of spatial infinity, disjoint from temporal infinity, is given by 
exponentiated versions of Eddington-Finkelstein coordinates. In our notation, this 
corresponds to 

lrr,l/2 -1/2 /Tnl/2 1/2 

where as usual lu denotes coordinates on S^. Here 

^hh,+ = {shh - = 0} 

is the characteristic surface given by /i = in T > (i.e. the front face of the blow 
up of the corner), and 

J^hh- = {sbh,+ = 0} 

is its negative time analogue. The change of coordinates (pbh,+, m) (sbh,+ , Sbh,-) 
is a diffeomorphism from (0,oo) x {0,6) onto its image, i.e. these coordinates are 
indeed compatible. As we show in the next section, the metric is C°° and non- 
degenerate on M, and the boundary faces Sbh,+ = and Sbh,- — are characteristic. 
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We can again extend M to M , which has only two boundary faces (the two temporal 
ones) by allowing Sbh,±, and analogously Sds,±, to take on negative values. Thus, 
M has six boundary faces, 

tf + , tf _ , .Fbh,± , -^MS.i: 

called the future and past temporal faces, and the future (+) and past (— ) black 
hole and de Sitter scattering faces. 




Figure 2. On the left, M is shown, while on the right its blow- 
down M. The time and space coordinate lines corresponding to 
the product decomposition are indicated by thin lines in the inte- 
rior. The temporal boundary hypersurfaces of M are continued by 
thin lines, as are the characteristic surfaces !Fh\x.± and Tds.±^ to 
show that the Lorentz metric extends smoothly across .?i5h,± and 
J-ds,± (but not across the temporal face!). The extended spaces 
are denoted by M and M. Valid coordinates near H .Fbh,- 

are (apart from the spherical coordinates) Sbh,+ and Sbh,-i as 
shown. 

The following propagation result follows directly from the properties of this blow- 
down. 

Proposition 1.2. If u satisfies Du = and has C°° Cauchy data on a space-like 
Cauchy surface E C M n {t > 0}, for example Y, — {t = 0} (i.e. Sbh.+ = Sbh.-A 
then w e C°°(M°). 

Combining Proposition 11.21 and Theorem 11.11 leads to the main result of this 
paper: 

Theorem 1.3. // u satisfies Du = and has C°° Cauchy data on a space-like 
Cauchy surface S C M H {i > 0} then there exists a constant c and e > such that 

(1.14) u-ceAl,^{M)^x^J\^,^{M) 

near the future temporal face, tf + . 
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Remark 1.4. Our methods extend further, for example to Cauchy data at t — 
which are conormal at dX, of growth m > —2, at the boundary. Standard 

hyperbohc propagation gives the same behavior at J^hh,+ and J^dS,+ in p < C, 
and then the resolvent estimates for the 'spatial Laplacian' Ax (described below) 
apply to yield the same asymptotic term but with convergence in the appropriate 
conormal space, including conormality with respect to J-^bh,+ and J-ds,+ - 

Remark 1.5. Dafermos and Rodnianski [4J have proved, by rather different methods, 
a similar result with an arbitrary logarithmic decay rate, i.e. an analogue of 

w-ce (logp)-^^Of^(M) 

for every N. In terms of our approach, such logarithmic convergence follows from 
polynomial bounds on the resolvent of Ax at the real axis, rather than in a strip for 
the analytic continuation; such estimates are much easier to obtain, as is explained 
below. 

As already indicated, by looking at the appropriate compactification, one only 
needs to study the asymptotics near tf+ in M (or equivalently, M). We do this by 
taking the Mellin transform of the wave equation and using high-energy resolvent 
estimates for a 'Laplacian' Ax on X. A conjugated version of this operator is 
asymptotically hyperbolic, hence fits into the framework of Mazzeo and the first 
author [7j, which in particular shows the existence of an analytic continuation for 
the resolvent 

R{a) = {Ax - a^)~\ lmcr<0. 
Here we also need high-energy estimates for R{a). 

The operator Ax has been studied by the second author and Zworski in [9 , where 
it is shown (using the spherical symmetry to reduce to a one-dimensional problem 
and applying complex scaling) that the resolvent admits an analytic continuation, 
from the 'physical half plane', with only one pole, at 0, in Im a < e, for e sufhciently 
small. Bony and Hafner in [T] extend and refine this result to derive polynomial 
bounds on the cutoff resolvent, X-R(c)X; X G C^{X), as \a\ — > oo in the strip 
|Imcr| < e. This implies that, for initial data in C^{X), the local energy, i.e. the 
energy in a fixed compact set in space, decays to the energy corresponding to the 
0-resonance. In our terminology this amounts to studying the behavior of the 
solution near a compact subset of the interior of tf+. Our extension of their result 
is both to allow more general initial data, not necessarily of compact support, and 
to study the asymptotics uniformly up to the boundary at temporal infinity. This 
requires resolvent estimates on slightly weighted L^-spaces, which were obtained by 
the authors in [5] together with the use of the geometric compactification M (or 
M). For this to succeed, it is essential that the resolvent only be applied to 'errors' 
which intersect dM in the interior of JFbh.+ and J-ds,+- This turns out to be a major 
gain since the analytic continuation of the resolvent (even arbitrarily close to the 
real axis) cannot be applied directly to the initial data. Thus essential use is made 
of the fact that once the solution has been propagated to the scattering faces, the 
error terms have more decay. 

It is then relatively clear, as remarked above, that if one only knew polyno- 
mial growth estimates for the limiting resolvent at the real axis (rather than in 
a strip), one could still obtain the same asymptotics, but with error that is only 
super-logarithmically decaying. This observation may be of use in other settings 
where such polynomial bounds are relatively easy to obtain from estimates for the 



WAVE SYMPTOTICS ON DE SITTER-SCHWARZSCHILD SPACE 7 

cutoff resolvent, as in [1], or analogous semiclassical propagation estimates at the 
trapped set, by pasting with well-known high energy resolvent estimates localized 
near infinity. This has been studied particularly by Cardoso and Vodev [3], using 
the method of Bruneau and Petkov 2, Section 3]. 

This paper is structured as follows. In Section [2] both the compactifications and 
the underlying geometry are discussed in more detail. The 'spatial Laplacian' and 
relevant resolvent estimates are recalled in Section [3] and in Section |4] the main 
result is proved using the McUin transform. 



2. Geometry 

In this section the various compactifications of de Sitter-Schwarzschild space are 
studied after an initial examination of the simpler case of de Sitter space. 

2.1. De Sitter space. We start with the extreme case of de Sitter space, corre- 
sponding to m = in (|1.1[) and (|1.2[) . to see what the 'correct' compactification of 
M should be. However, rather than starting from the static model, consider this 
as a Lorentzian symmetric space. De Sitter space is given by the hyperboloid 

zl + ... + zl ^ zl+^ + 1 in M"+i 

equipped with the pull-back of the Minkowski metric 

Introducing polar coordinates {R, 0) in (zi, . . . , z„), so 

R = yJzl + ... + zl = yrr^i^, e = R-^{zi, z„) e T = z„+i, 

the hyperboloid can be identified with Rr x Sg ^^ with the Lorentzian metric 



1 



(r2 + l)d0^ 



where d9^ is the standard Riemannian metric on the sphere. For r > 1, set x = t 
so the metric becomes 

(l+a;2)-idx2-(l + a;2)de'2 



An analogous formula holds for r < — 1, so compactifying the real line to an interval 
[0, 1]t, with T = a; = r-i for a; < i (i.e. r > 4), say, and T = 1 - |t|-\ t < -4, 

gives a compactification, M, of de Sitter space on which the metric is conformal to 
a non-degenerate Lorentz metric. There is natural generalization, to asymptotically 
de Sitter-like spaces M, which are diffeomorphic to compactifications [0, 1]t x Y 
of M,- X Y , where y is a compact manifold without boundary, and M is equipped 
with a Lorentz metric on its interior which is conformal to a Lorentz metric smooth 
up to the boundary. These space-times are Lorentzian analogues of the much- 
studied conformally compact (Riemannian) spaces. On this class of space-times the 
solutions of the Klein-Gordon equation were analyzed by the third author in [lO] . 
and were shown to have simple asymptotics analogous to those for eigenfunctions 
on conformally compact manifolds. 
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Theorem. Theorem 1.1.]; Set s±{X) = n^±-^ - X. Ifs+{\ys_{\) ^ 
N, any solution u of the Cauchy problem for □ — A with C°° initial data at t = is 
of the form 

u = x'+'^^K+ + x'-'^^K^, v±eC°°iM). 

If s^{X) — s_(A) is an integer, the same conclusion holds ifv^ e C°°(M) is replaced 
by = C°°(M) + x^+W-^- W \ogxC°°{M). 




Figure 3. On the left, the compactificatioii of de Sitter space with 
the backward hght cone from q+ = (1,0,0,0) and forward hght 
cone from = (—1,0,0,0) are shown. 17+, resp. ri_, denotes the 
intersection of these hght cones with t > 0, resp. t < 0. On the 
right, the blow up of de Sitter space at q+ is shown. The interior 
of the light cone inside the front face Sq^ can be identified with 
the spatial part of the static model of de Sitter space. The spatial 
and temporal coordinate lines for the static model are also shown. 

The simple structure of the de Sitter metric (and to some extent the asymptot- 
ically de Sitter-like metrics) can be hidden by blowing up certain submanifolds of 
M. In particular, the static model of dc Sitter space arises by singling out a point 
on§g-\ e.g. (Jo = (1,0,...,0) e §"-i C M". Note that {02,..., On) e M""^ are 
local coordinates on S"~^ near go- Now consider the intersection of the backward 
light cone from go considered as a point q^ at future infinity, i.e. where T = 0, and 
the forward light cone from go considered as a point g_ at past infinity, i.e. where 
T = 1. These intersect the equator T = 1/2 (here r = 0) in the same set, and 
together form a 'diamond', i7, with a conic singularity at g+ and g_. Explicitly f2 
is given by z| -I- . . . -f < 1 inside the hyperboloid. If q+, g_ are blown up, as 
well as the corner dil n {t = 0}, i.e. where the light cones intersect r = in il, 
we obtain a manifold M, which can be blown down to (i.e. is a blow up of) the 
space-time product [0,1] x B"-i, with B"-i = {Z G R"-i : jZj < 1} on which 
the Lorentz metric has a time-translation invariant warped product form. Namely, 
first considering the interior f2 of we introduce the global (in fl) standard static 
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coordinates [t, Z), given by (with the expressions involving x valid near T = 0) 

nn — 1 



3 Z = {Z2, ■ ■ ■ ,Zn) = X ^\/l+x'^{62,...,0n): 
Sinht = _j!l±^ = - (1 + ^2)(^2 ^ ^ ^ ^ ^ ^2))-l/2^ 



It is convenient to rewrite these as well in terms of polar coordinates in Z (valid 
away from Z = 0): 



r = Jz| + . . . + z2 = Jl + zl,. -zl=x + Ol + ... + 91, 



Sinht = ^^"+^^ = {x' - (1 + x2)(02 ^ ^ ^ ^ ^ ^2))-l/2 ^ _ ^2)-l/2^ 

Y ^1 ~ ^n+l 

- r-i(z2, . . . , z„) = (0^ + . . . + 0l)-^^^{92, 0n) e 
In these coordinates the metric becomes 

which is a special case of the de Sitter-Schwarzschild metrics with m — and A = 3. 

Lemma 2.1. The lift ofQ to the blow up [M;q-|_,(7_] is aC°° manifold with corners, 
n. Moreover, [fl, dflH {t — 0}] is naturally diffeomorphic to the C°° manifold with 
corners obtained from [[0, 1] x B""^; {0} x 9B"~^; {1} x 9B"^^] by adding the square 
root of the defining function of the lift of {0} x B"-i and {1} x B"-i to the C°° 
structure. 

Remark 2.2. This lemma states that from the stationary point of view, the 'right' 
compactification near the top face arises by blowing up the corner dM. x 9B"~i, 
although the resulting space is actually more complicated than needed, since the 
original boundary hypersurfaces of the stationary space can be blown down to 
obtain a subset of M. 

The fact that this approach gives p = as the defining function of the temporal 
face, rather than x (hence necessitating adding the square root of p to the smooth 
structure), corresponds to the fact that, in the sense of Guillarmou [5], M is actually 
even. 

Proof. Coordinates on the blow up of M near the lift of q+ are given by 



X, Qj = 9j, j = 2, . . . , n; 

these are all bounded in the region of validity of the coordinates, with the light 
cone given by ^2^=2 ®| ^ ^ (which is why the factor \/l + x^ was introduced), so 
the lift f2 of f2 is X]j=2 ®i — ^- Sj=2 ® j ~ ^ ^ non- vanishing differential 
where it vanishes, this shows that f2 is a C°° manifold with corners. Near the light 
cone, dfl one can introduce polar coordinates in 0, and use 

N 

X, r = (^ 82)1/2^ (t^2,...,c^„) = r-i(e2,...,e„)es"-' 

i=2 
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as local coordinates. On the other hand, blowing up the corner of [0, 1]t x B"~i, 
where T = e~^* for t> A, say, which is equivalent to (sinhf)"^ there, gives coordi- 
nates near the lift of T = : 

r, (sinht)^^/(l - r^) x^, uj. 

Thus, one almost has a diffeomorphism between the two coordinate charts, hence 
locally between the manifolds, except that in the blow up of [0,1]t x B"~i, the 
defining function of the temporal face is the square of the defining function of the 
temporal face arising from the blow up of M. This is remedied by adding the square 
root of the defining function of the lift of {0} x B"^^ and {1} x B""^, i.e. of x^, to 
the smooth structure, thus proving the lemma. □ 

It is worthwhile comparing the de Sitter space wave asymptotics, which is 

u = x''-^v+ + v+e C°°(M), e C°°{M) + a;"^i(loga;)C°°(M), 

with our main result. The fact that the coefficients in the de Sitter expansion are 
C°° on M means that on M, the leading terms are constant. Thus, the de Sitter 
result implies (and is much stronger than) the statement that u decays to a constant 
on M at an exponential rate. 

2.2. Blow-up of the space-time product. We now turn to the compactification 
of de Sitter-Schwarzschild space. It turns out that while this cannot be embedded 
into a space as simple as M the final setting is not much more complicated. In terms 
of de Sitter space, the difference is that while spatial infinity in [0, 1]t x B-^, blown 
up at the corner, can always be blown down, the same is not true for temporal 
infinity. 

In fact, the 'black hole end' r = rbh resembles the 'de Sitter end' quite closely, 
which motivates the construction in the de Sitter-Schwarzschild setting. There is 
a simpler construction, depending on the choice of a constant A > 0, which does 
not quite work because of some incompatibility between the two ends (which whilst 
very similar, are not quantitatively the same). With fi as in (|1.2p . compactify M 
by compactifying M into [0, 1]t as in (|1.5p - p.7p . so 

T+ = Ta.+ = in t > c, 

and compactify (rbh,?'ds) as [rbhjJ'dslj to obtain 

[0,1]t X [rbh,r-ds] x = [0,1]^ x X. 

Then blowing up the corners 

{0} X {rbh} X S^, {0} X {rds} x 

(and analogously at T = 1, i.e. T_ = 0), gives a space denoted M. Thus, a neigh- 
borhood U = Ux,+ of the 'temporal face' T+ = is diffeomorphic to 

(2.1) [0, e)p X [rbh, r^s] x p = Tx,+/^i. 

In the interior of the temporal face, where /i > 0, this is in turn diffeomorphic to 
an open subset of 

(2.2) [0,e)Tx(rbh,rds)x§2. 
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If the same construction is performed, but using the smooth structure on the 
the compactification of X given by a = M^^^, i-e. ^1/2, and 

fx,+ = tI^I = e-^\ t > c, 
then a neighborhood U of tf+ as above is diffeomorphic to 

(2.3) [0, e)p X [rth, rds]i/2 x §5, p = fx^+/a, 

where [^bh, ''ds]i/2 denotes that a has been added to the smooth structure (or 
equivalently (r — rbh)"^^^ and (rds — J")^^^ have been added to the smooth structure). 
The distinction between ()2.ip and (|2.3|) is the same as between [0,1] t x X and 
[0, l]f X Xi/2 (where T is defined analogously to T), namely the square roots of the 
defining functions of all boundary hypersurfaces have been added to the smooth 
structure. We denote the resulting space by Mi/2- 

The subtlety is that the de Sitter and black hole ends need different values of A. 
So what we actually need is to paste together U\^^^^^+^\{r < and U\^g^+r\{r > ri} 
for some ri , r2 , rbh < ri < r2 < r^s where Abh and Ads are chosen in a way that 
reflects the local geometry neat the two ends. The transition map in the overlap, 
where re (?'i, ^2), is given by 

(pbh.+ ,r,c^) ^ (pdS,+,r,L^), PdS,+ = /i^^^/^^"-Vbh!f 

In the overlap, where ^ Oso fi >—> y^-^dsZ-'^bh-i jg smooth, this is a polyhomogeneous- 
diffeomorphism, in the sense that it is polyhomogeneous in the local defining func- 
tions of the lift of T = (namely Pbh.+ and Pds,+). In particular, the front faces 
of the blow-ups have well defined boundary defining functions, namely fi, up to 
multiplying by a C°° non-vanishing function, so we consider the resulting space a 
polyhomogeneous manifold with corners, where the 'polyhomogeneous' (as opposed 
to C°°) faces are only the future and past temporal faces, tf+ and tf_. There is 
also an analogous construction for Mi/2- 

As indicated already, we also want a preferred defining function (up to taking 
positive multiples) x of tf-(- in order to measure the rate of decay at the temporal 
faces; this should be polyhomogoneous-equivalent to the local defining functions 
Pbh,+ and Pds,+. We take this to be of the form x — f{r)e~*, / > smooth for 
r G (7'bh,''ds). Comparison with (ll.9p - (|1.10p shows that we need to take /(r) = 
^-i/(2Abh) fQj. J. near rbh, fir) ~ ^i^i/C^^ds) fQj. near rds, hence x = pH^'^^'^'^ near 

the black hole boundary of tf+, and x — p^/^'^''^^ near the de Sitter boundary. 
Then, in particular, a neighborhood of tf+ is polyhomogeneous diffeomorphic to 

(2.4) [0,e), X [rbh,rds] x§2. 

We still need to determine the values A at the two ends. Writing T = T+ = Ta.+, 
the dual metric (which is the principal symbol of the wave operator, □) has the 
form 

G = Aa-^\^T^dl - a^dl - r-'^dl 

in the original product compactification, with dr = = 2[3d^. The change of 

variables from r to /i is smooth and non-degenerate, i.e. 2/3 = d/x/dr 7^ for fi close 
to 0, i.e. r close to rbh or rds. Note that /3 > for r near rbh, < for r near rds 
since ^ > for r e (rbh, ''ds). 

After blow-up, in the coordinates (p, /i,a;), 

G = 4p-iAV'9p' - 4M/32(ap - p-VSp)" - r-^dl 
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Thus, 

G = {AX^p'dl - m^^d, - pd,r) - r-^dl 

If we set A — /9(r'bh) > or A = — /3(rds) > therQ the p^9p terms cancel, so locally 
near rbh 

where 7 is C°° by Taylor's theorem, and there is a similar expansion at rds- Thus 
the choice of A determines the compactification M, and it is only at this point that 
the compactification has been specified. Note that this metric is a C°° Lorentzian 
b-metric on [0, e)p x (rbh — e, '"bh + e)r x §^ (i-e. is non-degenerate as a quadratic 
form on the b-cotangent bundle), in particular it is C°° across /i = 0. Denoting this 
extension of M by M (which is now non-compact); g becomes a polyhomogeneous 
conormal Lorentz metric on M, smooth near /i = (where there is a well-defined 
smooth structure). We write T for the set given by /i = 0, i.e. the boundary 
hypersurface of M that is no longer a boundary hypersurface of M . 

For this metric J- is characteristic, and one has the standard propagation of 
singularities in p > 0. In particular, for C°° initial data the solution is smooth 
in p > across = 0. In fact, writing covectors as + C,d^ + Yl,'nj di^j^ i-^- 
{p, fj,,uj,£,,CjV) si-re coordinates on ^T*M, the dual metric, considered as a function 
on ^T*M, is 

so the Hamilton vector field of G is 
Hg =8(7$ + f3'C)pdp - 8/32 (K - 

with H(^^ jj^ denoting the Hamilton vector field of the standard metric on the sphere. 
The conormal bundle N*{pi = 0} is /i = 0, ^ = 0, 77 = 0, so at this set 

and so is indeed tangent to N*{p, = 0}, and it is non-radial off the zero section 
(where C 7^ 0) as long as p 7^ 0. 

At dJ-, i.e. at p = 0, however there are radial points over the conormal bundle of 
J^. Rather than dealing with them directly, which can be done in the spirit of ^Gj, 
we reduce the problem to the study of the high energy behavior of the resolvent of 
the spatial Laplacian (which gives more, in fact), which was performed in 

2.3. Blov^r down of spatial infinity. We now discuss the manifold with corners 
M, in which spatial infinity is blown down. A valid coordinate system near the 
image of the black hole end of spatial infinity, disjoint from temporal infinity, is 
given by 

,^1/2 -1/2 /Tnl/2 Tnl/2 1/2 

where as usual lu denotes coordinates on S^. In these coordinates the dual metric 
is 

G = l{sl^,+dl^,^ + sl^^_dl^^ _)- 2(/3(rbh)2 + /52)9.,,„^a,,,,_ - r-^d^ 



^In terms of the 'spatial Laplacian', Ax, described in Section |3] — /3(r'bh)^ and — /3(r;js)^ are 
the asymptotic curvatures, hence they are natural quantities from that point of view as well. 
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and the boundary faces Sbh,+ — and Sbh.- = are characteristic. We can also 
extend G to a smooth non-degenerate Lorentz metric on M°. Recall that M contains 
M as a closed domain with corners, namely locally in M we simply allow the four 
boundary defining functions Sbh,±; Sds,±, to assume negative values (but we do 
not extend M across the temporal faces). This calculation shows the following 
important fact: 

Lemma 2.3. The d'Alembertian □ e Diff^(M°), and indeed □ e Difr6(M). More- 
over the scattering surfaces Thh,± = {sbh,^ = 0} and J-'ds,± = {sdS,T = 0} are 
characteristic. 

Proposition [L2l is an immediate corollary of this lemma and standard hyperbolic 
propagation. 

3. Resolvent estimates 

Next consider the 'spatial Laplacian', resolvent estimates for which constitute 
one of the key ingredients in our analysis. From ()1.4|) with T = e^*, it follows that 

□ {{TDrf - a^r-'^Dra'^r^Dr - a^r-'^A^) . 

Recall that not precisely T, but rather T^^^^ and T'^''^ , were used above to construct 
the compactification. 

By definition the spatial 'Laplacian' is 

Ax = a^r-'^DrO^r'^Dr + aV^^A^ 

Near a = 0, where a can be used as a valid coordinate in place of r, 

Ax = (3r-^aD^(3r^aD^ + a^r'^A^ e Diffg(Xi/2). 

This is not the Laplacian of a Riemannian metric on X; however it is very similar 
to one. It is of the form d*d with respect to the inner product on one-forms given 
by the fiber inner product with respect to the 'spatial part' H = o?d1 + r~^9^ of 
G and density on X given by dh — a~^r^ dr dto. In what follows we will also view 
Ax is a 0-operator on self-adjoint on 

L^{X, \dh\), \dh\ = a-^r"^ \dr\ \duj\ = a-^\f3\~^r^ \da\ \dLu\, 

and we will use the techniques of [7] to study its resolvent. It is also useful to 
introduce the operator 

L — aAxa~^ , 

which is self-adjoint on 

L^{X,a-'^ \dh\) = aL^{X,dh), a-^\dh\ = a-^j^j-V^ \da\ \duj\. 

Thus, this space is Lq{X) as a Banach space, up to equivalence of norms. 

Let a = a^/'*'*''' S C°°{X), a > 0, for r near rbh or a = 0;^/'*'''^ for r near r^s- The 
normal operators No±h{L), No^dsiL) of L in DiffQ(Xi/2) at r = rbh, resp. r = rds, 
are 

NoML) = Agh%bh(AH3), iVo,ds(i) = A^s^o,ds(AH3), 
where Aga is the hyperbolic Laplacian, explaining the usefulness of this conjugation. 
In particular, it follows immediately from 7 (with improvements from [5 ) that the 
resolvent 

7^(c^) = {L- a^y^, on Lg(Xi/2), Ima < 0, 



14 



RICHARD MELROSE, ANTONIO SA BARRETO, AND ANDRAS VASY 



continues meromorphically to a strip | Imcrl < e as an operator between weighted 
L^-spaces (as well as other spaces): 

7^(a) : a'Ll{Xy^) ^ a-'Ll{Xy^), 5 > e; 

we keep denoting the analytic continuation by TZ{ij). Thus, 

R{(t) = (Ax - 0-2)-! = a-^n{(j)a on L^{X, \dh\), Imcr < 0, 

continues meromorphically to a strip | Imaj < e 

R{<t) : a^L^iX, \dh\) a^^L^iX, \dh\), 5>e. 

The result we need is proved in [8j, giving polynomial bounds on the resolvent 
in a strip around the real axis. 

Proposition 3.1. If e > is sufficiently small the only pole of the analytic con- 
tinuation of the resolvent R{<t) in Imcr < e is a — Q, which is simple, with residue 
given by a constant 7 and for each k and S > e there exist m > and C > and 
M such that 

(3-1) \\a^'"^R{(T)\\c(a-^&'>H^{x^,2),c'-{x)) < C\a\^\ 

for \a\ > 1, Imcr < e. 

4. ASYMPTOTICS FOR SOLUTIONS OF THE WAVE EQUATION 

We now proceed to study the asymptotics of solutions of the wave equation at 
tf+. 

Suppose It is a solution of the wave equation, Du — 0, and u is smooth on 
M°. Energy estimates show that u is necessarily tempered, in the sense that u E 
p~^L^{M) for some s > 0. This is shown directly below. 

Let (j) = 0o(p) e C°^{M) be a cutoff function, with p e C^{[0, 1)), identically 1 
near 0. If u satisfies Du — and u is smooth in /i (i.e. across the side faces), then 
(pu is smooth in /i, and 

□ (0u) = [□, cl,]u = f 

where / is also smooth in /z, and vanishes in a neighborhood of the temporal face 
- since [□,</>] € DiffJ,(M), by Lemma [2.3[ is supported away from the temporal 
face. Moreover, v = (pu is the unique solution of Dv = / in M° with v = for p 
sufhciently large. 

Now we wish to take the Mellin transform in T = e^*, for functions supported in 
a neighborhood of the temporal face, namely in U = {p < 1}. Such a neighborhood 
is equipped with a fibration AI — ^ X, extending the fibration Af X in the 
interior, and there is a natural density \dt\ = on the fibers. In coordinates 
{p,a,u}) = {T^^^ /a,a,Ljj) valid near the temporal face boundary at r = r^h this 
density takes the form So, the Mellin transform can be taken with respect to 

this fibration and density. Thus the map v v from functions supported near the 
temporal face to functions on x X , fl C C, 

v{a,z)^ j r%(r,z)^. 

If V is polynomially bounded in T, supported in T > 0, with values in a func- 
tion space 7i in z, this transform gives an analytic function in a lower half plane 
(depending on the order of growth of v) with values in Ti. 
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In fact, writing the integral in coordinates vahd near the boundary of the tem- 
poral face, p — T'^^^^/a, 

The integral is then the Mellin transform of v with respect to p evaluated at tr/Abh- 
Thus, if V is smooth on A/°, supported in {p G I}, I C (0, 1) compact, i.e. near but 
not at the temporal face, then it is in fact analytic in C with values in functions of 
the form a'''^C°°, with C°° seminorms all bounded by Ck{a')~^ , k arbitrary. If v is 
just supported in p < 1 and is conormal on M, then v is analytic in a lower half 
plane with values in functions of the form a^'^C°° . 

Assuming for the moment that <j)u is polynomially bounded in T, = / 

becomes Nb{U\)(f)u — f, where 

Nb{D) = a-^{a^ - Ax) = a'^ {a^ - a^r'^ Dra^r^ Dr - a^r'^A^) , 

so 

(a2 - Ax) 0u = 

If (f>u is polynomially bounded in T, then both / and (j)u are analytic in Im cr < 
— C, and as / is compactly supported in p, f is in entire analytic with values 
in functions of the form ff'^C°°, with C°° seminorms all bounded by Ck{(j)~^ , k 
arbitrary. Thus, 

(pu = R{a){a^f), Im cr < -C, 

and we recover (j)u by taking the inverse Mellin transform. 

We now return to arbitrary (not a priori polynomially bounded) u, f — □(0u), 
as above. Thus, / is analytic in all of C, with values in functions of the form 
a^°'C°°, with C°° seminorms all bounded by Ck{(j)"'^, k arbitrary. Now note that 
the inclusion 

is continuous for every s > 0. Thus, the inclusion 
is continuous if 

A(e + 6) <2, A = max(Abh, Ads), 
and Imcr < e, which is to say 

a'-a^C°^{X,/2) a-^a'H^iX,^^) 

is continuous. 

In particular, then 

for all k in Imcr < e < S (with new constants), < e < 6 sufficiently small. 
Proposition 13.11 shows that, for e > sufficiently small and for all N and k, 

(4.1) l|a"'"^M(a'/)llc«m < Ck\a\-^lma < e. 

The inverse Mellin transform of w = R{a){a^f) is 

w{T,z) = {2n)-^ [ T-"'w{(j,z)da. 
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Thus, 

In view of (j4.ip (in particular the analyticity of w in the lower half plane with the 
stated estimates), w = for T < (as can be seen directly by shifting the contour 
to Ima — — C, using the off spectrum resolvent estimate \\R{<j)\\c(L2(^x,\dh\)) < 
I Im(cr^)|~^ and letting C —> +oo). Since the unique solution of □« = f ,v supported 
in T > 0, is it follows that w — cjju. 

Shifting the contour for the inverse Mellin transform for w to Im cr = e gives a 
residue term at 0, and shows that 

p-^{cj,u-v)eL^{[0,S)p;C"°{X)), 

where v arises from the residue at 0, hence is a constant function. Note that (at the 
cost of changing e) this is equivalent to the analogous statement with p replaced 
by p. The derivatives with respect t satisfy similar estimates. Hence, the same 
estimates hold for the conormal derivatives with respect to p (or equivalently p). 
We thus deduce the leading part of the asymptotics of u at the future temporal 
face, tf+. 

This completes the proof of Theorem ll.il The main result. Theorem II .31 follows 
from the combination of this result with Proposition 1 1.21 
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